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Abstract 

We predict several effects associated with the optical response of systems prepared in a nonequi¬ 
librium state by impulsive optical excitations. The linear response depends on the phase of the 
electric field even if the initial nonequilirbium state has only populations, no coherences. Ini¬ 
tial quantum coherences induce additional phase dependence which also shows new resonances in 
nonlinear wave mixing. In systems strongly driven by an external optical field, the field frequency 
generates a phase dependent probe absorption. This gives further control to manipulate the relative 
contribution to the linear signal due to initial populations and coherences. 
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I. INTRODUCTION 


Spectroscopic signals are commonly expressed in terms of the po 


the molecule by external driving 


a 


arization induced in 
, 3], the polarization 


For n -wave mixing processes 
is calculated perturbatively to n — 1-th order in the field-matter interaction and written as 
a combination of n point matter correlation functions which depend on the initial state 
of the molecule. In the frequency domain, the polarization is expressed in terms of the 
n-tli order optical susceptibility y'"F For systems initially at equilibrium, due to time- 
translational invariance symmetry, x (n> only depends on n independent frequency variables. 
This is because the sum of all n + 1 field frequencies must add up to zero. This reflects the 
energy conservation of the field. 

If the system can be prepared in an arbitrary nonequilibrium state the time-translational 
invariance symmetry is broken and the nonlinear optical response is governed by a more 
general y- n) that depends on n+1 independent frequency variables. This one extra frequency 
variable yields new resonances that are not present in the n-tli order equilibrium response. 
To make a consistent connection with n-wave mixing results obtained with equilibrium state 
preparation, we introduce a generalized n-th order optical susceptibility x (n ' 1 , which depends 
on the nonequilibrium initial state. This initial state can contain populations and coherences. 


Such nonequilibrium state prepared 
femtosecond Raman spectroscopy 4 


yy an ultrashort actinic pulse is commonly studied in 
7]. In a different context, systems prepared in the 


nonequilibrium superposition of quantum states induced by strong monochromatic laser 


field have been extensively studied in quantum optics B 0 and quantum thermodynamics 
10-12]. Monitoring field induced quantum coherence by the transmission of a weak probe 


shows many interesting effects including electromagnetic induced transparency 


i 


n 


Q, 


lasing 


without population inversion [1J, LUa] and coherent population trapping [mini. Molecules 
in open junctions is another example of nonequilibrium preparation, where the molecule 
subjected to a chemical or thermal bias, across its boundaries, can be in a stationary state 
10] For such initial conditions the standard fluctuation-dissipation relations are violated 
18-20]. In all the above examples nonequilibrium preparation affects the response of the 
system with respect to an optical perturbation. Typically the optical response of a system is 
probed by nonlinear n -wave mixing. For instance Raman y® measurements jfi] can detect 
phase and amplitude information about resonances. In this paper we show that, similar 
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information about phase can be obtained from simpler linear measurements described by 
yd) w hich depends on two independent frequencies. The linear signal now depends on the 
phase of the electric field even if the system is initially prepared in a population state. 
The initial coherence further results in additional phase dependence which can be used in 


coherent control schemes 


21 


I 22]. Another consequence of the nonequilibrium preparation 
are the new resonances in n wave mixing. These are at the frequencies shifted from the usual 
n -wave mixing phase matching condition by an amount determined by the phase of the non 
stationary state corresponding to initial coherence { 23 ]. We further study an example where 
the initial conditions are considered explicitly for a three-level model systems driven by a 
strong monochromatic light source. We show how the control parameters of the driving 
field affect the phase dependence of the linear signal by manipulating initial coherences and 
populations. 

The paper is organized as follows: In section II, we calculate linear absorption of systems 
initially prepared in an arbitrary superposition state and discuss the phase dependence due 
to the field envelope. Compact expressions are derived by employing superoperator notation. 
In section III, we extend the calculation to nonlinear wave mixing signals and examine new 
resonances that show up because of the initial quantum coherence. In Section IV, we repeat 
the calculation for linear absorption of a strongly driven three level model system. The phase 
dependence now appears in the signal through the field frequency for both initial population 
and coherences states. We further examine the role of driving field parameters to control the 
relative contribution to the signal coming from population and coherences . In Section V, 
we summarize our findings and discuss possible experimental conditions where these effects 
can be realized. 


II. PHASE DEPENDENCE OF LINEAR OPTICAL SIGNAL 

We consider a multilevel quantum system driven by a classical optical field and is de¬ 
scribed by the Hamiltonian 

H(t) = H 0 + H int (t), (1) 

where H 0 is the system Hamiltonian and H int (t ) is the light-matter interaction 

tfmt(i) = V£(t) = V(€(t) + £*(t)). (2) 
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Here V — p + /T is the dipole operator which is partitioned into lowering (//) and raising 
(/A) operators and is responsible for the de-excitation and excitation between the molecular 
states, respectively. £{t) and £*{t) are positive and negative frequency components of the 
total electric field £(t) = £(t) + £*{t). We will calculate signals for an arbitrary initial 
molecular state. We employ superoperator notation where with each Hilbert space operator 
A, we associate two superoperators, denoted as A L (left) and A R (right), defined through 
their action on Hilbert space operator X as AlX = AX, ArX = XA. We further define the 
linear combinations A + = (A L + A R )/2 and A_ = A L — A R . A +(—) operation in Liouville 
space corresponds to an anticommutation (commutation) operation in Hilber t sp ace. Using 
this notation the frequency-dispersed heterodyne detected signal is given by |2J, [25] 


p r f°° / 

%r) = -Im[r(w) J dte i ^ t - t ^{TV L (t)exp( 


^ J Tt 


TO 


P(TO)J 


where T represents the set of parameters of the incoming fields. These will be specified 
later. The angular bracket (• • •) represents an average over the initial density matrix of the 
molecule prepared at time To- The external fields have finite envelopes centered at to > tq. 
The superoperators in Eq. (ED are in the interaction picture with respect to the free system 
Hamiltonian H 0 i.e., for any operator A, A v [t) = exp (iH 0 _(t — To))A u exp(—iH 0 _(t — To)), 
v = L,R. We further define the retarded Liouville space evolution operator Q(t — t 0 ) = 
(— i/h)6(t — To) exp [— %Ho-(t — To)] and the advanced evolution operator Q\t — ro) = 
(■ i/h)6(r 0 — t ) exp [—— To)]. In the frequency domain these propagators are written as 
G(u) = dte^-^Git - T 0 ) = \(wl - \Ho- + ie)- 1 and &(cu) = dte^-^Q^t - 
To) = \{wl — \H 0 - — ie) -1 . Here I is the Identity operator in Liouville space. T is the 
time-ordering superoperator 


TA^t^Bv'fo) = 6»(U - t 2 )A v (t 1 )B l/ >(t 2 ) + 0(t 2 - t^B^^A^B), u, v' = L , R. (4) 


The linear signal is obtained by expanding the exponential in Eq. (ED to first order in the 
field-matter interaction 


S^\uj-t 0 ,T 0 ) = -Im 
n 


r roo rt 

ih dt dne^-^ £*(oj)£( T i ) (V L g(t - r^V-Gin 
J— oo J— oo 



( 5 ) 

Note that To) = —jJ if p(to) is a stationary distribution (equilibrium or steady state). 

Otherwise it describes the evolution of the nonstationary initial state. The presence of 
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6{r\ — To) in Q{r\ — To) allows to extend the lower limit of integration for t\ to —oo. Making 
the change of time variable t — T\ = t\ , we write 


dtie lu ^ to) £*(u)£(t-ti-A,o){V L G(ti)V-G(t-ti-To)) 


p( t o). 


s {1) (u-t 0 , To) = — Im | ih J dt 

( 6 ) 

Performing the inverse Fourier transformation for G(t—ti—To), £(t—ti —to) and integrating 
over ti and t, we obtain 


S^\w, to, To) = - Im 
n 


ih 


dcjQ 

2n 


rf, / 

—£-(w)£M)(Vi SK + u.„)V_ 6 (wo)} p(to) 


2 t 


x 27 T( 5 (o; — ojo — uj[) 

The cu'o integration can be carried out easily which results in 


S ,( 1 ) (o;;to,To) = ^Im 


_1 e -,(„-„i)(.„- T o) . (8) 


with the system while the coherence is alive. 


photosynthetic complexes 


which gives 


f 

J —c 


d(jj[ 

2n 


S(u[)Gab{u - U>i)e 




= — — lim 


du>[ 


h >o ,/_ oo 27r 


dt£(t) — 


e *OJj(t-To+t 0 ) 


i 


(Uj - (W - W a6 ) - 777 
dt 0(* - To + to) £(*) eh^-^iCt-ro+to)^ 


( 9 ) 


S {1 \uj ,; to-To) = ^ Im 


ab 


did(t+t 0 -T 0 ) £%lj)£(^(V l G(cj)V-) pa e liu - Wab)i e- iu ^ to - To) 

( 10 ) 


( 7 ) 


This expression holds for an arbitrary initial density matrix. 

We now assume that the molecule is prepared in a superposition of eigenstates p(r 0 ) = 
^2 fl b P»i> I °b)) by a weak impulsive pulse centered at To and calculate the linear absorption 
of a second weak probe pulse. We neglect dephasing and assume that the probe interacts 


'his is the case in recent experiments in 


26] and solar cells [27]. In Section IV we consider a particular 


example where the strong driving held induces quantum coherence and it is long-lived. For 
a superposition state ab)) the contour integration over uj\ in Eq. jED can be carried out 


where in the second line 77 is an infinitesimal positive number used to satisfy the causality 
condition for the retarded propagator. In last line we use the definition of Heaviside theta 
function. The signal is finally given by 


5 





















In the following we discuss various limits of the signal and its dependence on the phase of 
the field. 


A. Systems at equilibrium 


The standard equilibrium result is recovered from Eq. m in the limit To —> — oo where 
the molecule initially may reside in some population state | aa)) (no coherence) or, in thermal 
equilibrium i.e., p eq ( r (J ) = e~P H ° /Z, where ft is the inverse temperature and Z = Trje - ^ 0 ] 
is the partition function. The linear signal then reduces to 



h 


Im 




( 11 ) 


where x^(— u>;oj) = (VlQ( co)V-) Peq which depends on a single frequency and £(u>) is the 
Fourier transformation of the total electric field £ (t ). Making the rotating wave approxima¬ 
tion (RWA) and assuming that the molecule initially is in a ground electronic state |a), we 
get for the linear signal 

'S'eq’M = 




|P 

ca | Paa 
uj — 4- in. 


( 12 ) 


This solely depends on the power spectrum of the field \£{uj)\ 2 and is independent of its 
phase. 


B. Nonequilibrium state and long pulses 


When the pulse envelopes are long (continuous wave (CW)) i.e., the field £(t) = £ e * u ’ 11 , 
and £{oj) = 2t£8{(jJ — uq) the integrated signal (Eq. (JTQD) is given as 


SV = 


duj s (1 \uj) = | J2 lm \ £ *£( v L g M v - 


2n 


ab 


' Pab 


Gabi^l ~ U>l) 


(13) 


The signal is now independent of t 0 and r 0 . Note that even if the system is initially prepared 
in a superposition state, i.e. there is nonzero initial coherence p 0 & the CW signal does not 
show uj a b dependent resonances. 


C. Initial nonequilibrium state with populations and coherenes 

Starting with a general nonequilibrium state with initial populations and coherences and 
following Eq. ca with the assumption that the field envelopes are centered at t 0 = r 0 we 
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obtain the linear signal 


S^\oj) = — Im[ 
h l 


hm 

n 


£ £»£(w - wj 

ab 


( 14 ) 


where 




<&£(*) e^, 


( 15 ) 


J 0 

is the one-sided Fourier transform of the electric field which makes the signal dependent on 
the phase of the field. We define the generalized linear susceptibility 


X {1) (-uj,uj[) = ^(V l Q(uj)V^) p^Slu-u'i-Uab), ( 16 ) 

ab 

which depends on two independent frequency variables. Expansion of matter correlation 
function in molecular eigenstates gives 


<y L g(u)V-) M 


^ ^ Pabpca 
c 


pbc 

-Uj — cj C 5 if] 


pcb 

M ~ ^ac + 


( 17 ) 


where /i ac = fi ac + /i* c is the matrix element of the total dipole operator. 

The generalized susceptibility in Eq. m may also be recast in terms of forward (Q) and 
backward (Q^) Liouville space propagators 

X ( 1 ) (-cn,n;() = ■£- ^ Pab((I\V L G(uj)V- [Q(u - uj() - (uj - uj[)\\ab)) . ( 18 ) 


In the to = h) limit we could also write the signal directly from Eq. (8) as 


5' ( - 1 ^(a;) = — Im 


ih 


duj[ 
27r 


£*(u)£(u[)(V L g(u)V-g(u - u[)) 


p( t o) 


( 19 ) 


This equation contains an additional Green’s function Q(uj—uj() compared to systems initially 
at equilibrium. This frequency-domain Green’s function is given by a one-sided Fourier 
transform of a corresponding time-domain Green’s function in Eq. (6), which depends on a 
superoperator time ordering via the Heaviside theta-function. The latter is an intrinsic part 
of the bookkeeping of the field-matter interaction which is typically done by assigning the 
time-ordering operator to the matter correlation function, while allowing the fields to evolve 
in an unrestricted manner. The frequency domain electric field £(cu() is given by a complete 
Fourier transform of the time-domain field envelope. Alternatively, the burden of time¬ 
ordering can be placed in the field correlation function while allowing matter to evolve in 
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time without restriction. In this case the complete Fourier transform of the Green’s function 
results in the delta-function as shown in Eq. fUHl) . while the electric field is transformed via 
one-sided Fourier transform as in Eq. (fHl) - (fT51h Both representations are identical and 
the choice can be made according to the convenience in bookkeeping of the field-matter 
interactions and details of the signal. When time-translational symmetry is restored for an 
initial equilibrium or steady state (with or without coherence) the Green’s function Q(cu — 
(ju[) = ) and the linear signal or the susceptibility yd) depend on to a single frequency 

variable. The breakdown of time-translation symmetry due to initial nonequilibrium state 
gives rise to signal described by %d) which depends on two frequency variables as shown in 
Eq. m- 


1. Phase dependence induced by initial populations 


From Eq. (Tl4l) we separate the contribution to the signal from initial populations 

2 
h 


= ilm[j]£*( W )£(u>){Ue(u.)r_> ta ]. (20) 


Despite the fact that initial state contains only diagonal elements of the density matrix 
(populations), the signal in Eq. fl20l) depends explicitly on the phase of the electric field. 
This is because the initial state at r 0 is not an equilibrium state which makes the signal 
depend on 8(co) instead of £(cu). It is only in the limit To —» —oo that the phase dependence 
disappears. 


To demonstrate this phase dependence we consider a linearly chirped [281, 29] Gaussian 


pulse with spectral phase <p(uj) = 0 O + ~ <w) 2 /2. The electric field is given as 


£(u) = V7r^oe“ (w “ s>e)ai “ /4 e^" (w “ ac)a/2 , 


( 21 ) 


where T 0 /\/2 is the field-transform-limited temporal width, 5u>l = \[2jTo is the spectral 
width, u) c is the central frequency of the field and <fi" is the quadratic phase inducing the 
linear chirp. The corresponding temporal profile is 

£(t) = (22) 

where p = xT — an< ^ V^o = Tq. The temporal width of the chirped pulse is given 
as T p = T 0 y/l + (20"/Tq) 2 and the instantaneous frequency is u>(t) = Q c + 2 at with a = 







(23) 


( 2 4>")/\Tq + ( 2 4>") 2 ]. Eq. (fTHD is then given by 


£(u) 



(u-u c ) 2 Tg/4 e i<t>"(u-u c ) 2 /2 M ^ g r g 


~UJ ~ Ld c ~ 

- 2VT - 
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where Erh [z] is imaginary error function. It is then clear that the chirp rate dependence in 
the signal enters through £*(oj)£(oj) oc Erh ^^h ] function. 


2. Phase dependence due to initial coherence 


From Eq. (1141) the contribution of initial coherences to the signal is 

sShM = \ lm -[ Y £*(u)£(u^uj a b)(v L g(uj)V-) Paib 

a : b : a^b 


(24) 


The signal now shows a new resonance which depends on the initial coherence frequency oj a i, 
between levels a ans b. The field envelope depends on the shifted frequency £{uj—uj a b) which 
also generate a phase dependent signal. For a chirped pulse, the chirp rate dependence in 
the signal is £*(u ))£{uj — uj a b) oc Erfi[ ^~^? ]. 


3. Application to a three level model system 


In the following we demonstrate the phase dependence in the linear signal for the three 
level model system shown in Fig. [D(a). We assume that the two lower states |a) and | b) are 
initially in a maximally coherent state i.e., p aa = pbb = p a b = Pba = \ and all other elements 
of the density matrix vanish. Dipole transitions are allowed between states |a) —» \c) and 
| b) —? |c). Using Eq. (HU) and making the RWA, the signal can be written as Sl'J( oj) = 
J2i, j=a ,b S ij\u) where 



— Im 
n 


£*(uj)E(tjj — ujij) 


UJ — UJ c j + XT]. 


(25) 


In Fig.[T](b) and Fig.[T](c) we display linear transmission vs the chirp rate and frequency when 
the system starts in population states p aa and pbb respectively. SaJ and show peaks 
corresponding to uj ca = 0.8 eV and uj c b = 0.7 eV transitions respectively. The population 
oscillations in the signal come from the error function which oscillates faster with the chirp 
rate <p" for higher oj—Cj c . We use the central frequency uj c = 0.5 eV which implies that for oj co 
transition the oscillation is higher compared to a : c b transition. Also the signal is asymmetric 
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FIG. 1. (a) Level scheme for a three level model system with two lower states | a) and \b) and one 
upper state | c). The dipole transition is allowed between the states | a) \c) and \b) |c). The 

states \a) and \b) are initially in a maximally coherent state i.e., p a a — Pbb — Pab — Pba — \ • (b) 2D 
plot for the linear transmission signal vs the detected frequency uo and chirp rate starting 
with initial population in state |a), (c) Plot for with initial population in state \b) and (d) 
the combined signal due to population Spop = (e) The signal S^ h = due 

to initial coherence ( pabipba )• (f) The total signal S^H including both population and coherences. 
The parameters are given as follows : Ub a = 0.1 eV, cc CCi — 0.8 eV, Oo c — 0.5 eV, p = 0.004 eV, 
T 0 = 6.6 fs. 


with respect to the chirp rate (j)". In fact, it is stronger for positive (/)" as compared to the 
negative one and is due to the amplitude of the field envelope which is higher for positive 
chirp. In Fig. [T] (d) we display the total signal due to initial populations Spop = Sm + Sj^. 
In Fig. [T| ( e) we plot the signal + Sj^ due to initial coherence. This signal appears 

to be stronger compared to that due to populations even for negative chirp. In Fig. [Iff) we 
display the total signal S^t due to initial population and coherence. 
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III. NEW RESONANCES IN NONLINEAR WAVE-MIXING 


The Above results for the linear signal can be easily extended for nonlinear signals. In 
the following we present expressions for three wave-mixing (TWM) and four-wave mixing 
(FWM) signals from an initial superposition state. The electric field 8 is a sum of two 
(for TWM) or three (for FWM) monochromatic field modes Ei,i = 1,2,3. Because all the 
frequency modes overlap in time, we must include all possible permutations of these modes 
to calculate the full signal. We first write the TWM signal in terms of the total field 8 


S<2) ^ = ! Im / Tk S 


(26) 


where the second-order generalized susceptibility is 




(27) 


ab 


This can be alternatively recast as 




ab 


(28) 


Expanding the matter correlation function in system eigenstates gives 


(Vi Q(u)V.Q(u - u.;)Vi) p „, = 


cd 


l^dcl^caGcbi^ kfi) ( l^bdGdb{^) l^dbGcd{^) 


T l^da/J'cbGaci^ ^l) ( l^dcGadi^') l^cdQdci^^^) 


(29) 


Extending these results to FWM we can similarly write 
r00 dco[ f°° dco 2 rco ' h ■' 


S^\oj) = — Im 
h 


d(jJ O 


I —oo J —oo J —oo 


£\u)£(uj[)£(u' 2 )£(u' 3 )x {3 \-u,uj , 1 ,uj' 2 ,u , 3 ) , (30) 


where 


X^\-u,u[u' 2 ,u' z ) = ^{V L G{u)V-G(u-u[)V-G(u-u'i-u' 2 )V-) Pab Siu-u^ -oj 2 -u' z -u ab ), 

ab 

ab 

[£/(^— oj'i — co 2 —(Ug) — (cu—— co 2 — W 3 )] | ab)). (31) 
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Expanding in eigenstates finally gives 


{v L g(u)v.g(uj - u[)v-g(u- = 


5 Z pab 


l^cal^dcl^bel^edgcbi,^ kfi ^2) Odb(^ kfi) ^Seb{^} 


cde 


T l^cal^bdl^ecl^de Gad{^ ^1 
T l^ecl^def^bdl^ac Gcb{^ ^1 

“ 1 “ l^eal^del^cdl^bc Gacif^ ~ 


^2) Gcdi^J ^l) ^ Ged (^) ^ce(^)^ 

CJ2) Gcd{^ ^l) ^ 'Ged( 0 ce(^)^ 

— CJ2) Gadis fj J ^l) ^^ae(^) 


The equilibrium limit for TWM and FWM signals can be obtained in a similar way as was 
done in Eq. (HD- 

Below we present the FWM signal for the three-level model system, shown in Fig jlk). 
The complex field amplitude contains three monochromatic modes and a spectrally broad 
Gaussian probe field. £(t) = y ^ =1 £j exp(?'kp'j — iujjt) + f ^£ 4 (cj) exp(ikr — loot) where 
the probe field £ 4 ( 0 ;) = y^exp [(cn — cJ c ) 2 /2a 2 ]. We select the following phase matching 
direction for the signal k = ki—k 2 +k 3 . The signal is calculated using Eq. (1301) . As before we 
assume that the system is in a maximally coherent state. Various Liouville space pathways 
and corresponding expressions for the FWM signal are given in Appendix A. In Fig. [21(a) 
we display the signal starting with a population Spop(A') vs the detuning A = u — uj\ + 
u 2 — 0 J 3 . Various peaks correspond to the different quantum pathways of matter represented 
by frequency-domain propagators, which depend on some frequency combinations of the 
incoming field. The A = 0 peak corresponds to a detected frequency oj = 1.35 eV and is 
marked by an arrow. For the system initially prepared in a coherence the original peak at 
A = 0 splits into two peaks at A = Toy,/,, as marked by two arrows in Fig. [21(b). The spectra 
in (a) and (b) are not identical, as different matter pathways contribute to the signal if the 
system is initially prepared in a population or in coherence. In Fig. EKc) we compare the 
total signal S^(cj) = Sl%(cj) + S^l(cu) with the signal coming only due to the population 

/o\ 

iSpopfio')- The effect of initial coherence in the total signal is significant and many single 
photon resonances appear much stronger in the presence of coherence. In Fig. El(d) we 
compare the signal obtained for a system initially prepared in a non equilibrium population 
state vs equilibrium state. Note that the data for the population is scaled. The signal 
contain the same number of peaks but with different magnitude. The difference is solely 
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FIG. 2 . (Color online) Plot for four FWM signal vs the detuning from the phase matching A = 

(o\ /o\ /o\ /o\ 

l o —cJi +u )2 —^3 for different initial conditions. Sp 0 p - (a), S^- (b), *S t v 0 ^ al (black) and Spop(red) - 
(c), and Spop (black) and Seq (red) -(d). A single arrow in (a) and two arrows in (b) correspond to 
the splitting of a singe peak in population in two new resonant peaks in the presence of coherence 
as A = 0 is replaced by A = Any,. The parameters are given as follows: uj a = 0 eV, cj/, = 0.4 eV, 
uj c = 1.2 eV, ui\ = 1.1 eV, ui 2 = 0.75 eV, W 3 = 1.0 eV, = 0.5 eV, a = 10 eV, 17 = 0.002 eV. 


due to the non equilibrium preparation of the initial state. The two spectra coincide in the 
limit To —> — 00 whereby £ —> £. 


IV. LINEAR RESPONSE OF A STRONGLY DRIVEN SYSTEM 

A. Dressed state description for a strongly driven three-level system 

So far we did not specify how the system has been prepared in the stationary superposition 
of quantum states. We now consider a specific type of preparation using a strong driving 
field. We examine the linear absorption of a three level system, shown in Fig. (j3J), which 
is first prepared in a non-stationary state by driving the two lower levels |a) and | b) with 
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FIG. 3. Level scheme: Three level system with two lower states | a) and \b) and one upper state |c). 
The states | a) and \b) are driven by strong monochromatic field with frequency c^o- The optical 
transition is allowed between levels | a) \ c) and \b) |c). 


a strong monochromatic field of frequency ccq- The system is further in contact with a 
thermal bath that causes dephasing and relaxation. We consider a linear measurement 
taking this non-stationary state as an initial state. The system plus bath is described by 
the Hamiltonian 


H(t) = H 0 (t) + H B + H SB , 

H 0 (t ) = ^ hwi\i)(i\ - \e lUot \a){b\ + e~ lUot \b)(a\], 

i=a,b,c 

H b = ^2, 

k 

H S B = h ^2 al\i)(j\ + h.c (32) 

i^j=a,b,c,i<j 

where H^it) is the Hamiltonian for the system including the interaction due to the driving 
field. ak(a[) is the Bosonic annihilation (creation) operator for the bath, represented by the 
Hamiltonian H B and Hsb is the system-bath coupling Hamiltonian. The driving electric 
field £d(t) = £ 0 cos(u> 0 t ) and ujq is near resonant with the low energy states |a) and |6). The 
Schrodinger picture evolution of the reduced density matrix in the lab frame is given by the 
Bloch equations p{t) = C(t)p(t ) (see Appendix B). The driving field couples the population 
of the two low energy states i.e., p aa and pbb with its coherence p a b and maintains a finite 
long-time coherence. 

We shall probe this system prepared in a nonstationary state at time r 0 with a weak 
optical field £(t) which allows transitions between the states |a) —>■ \c) and \b) \c). We 
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FIG. 4. Ladder diagrams for a three level model system driven by a monochromatic field. Diagrams 
[(al)-(a8)] are for initial conditions starting with population states. Diagrams [(a9)-(al2)] are for 
initial conditions with coherence states. Time is increasing from bottom to top. 


write the light-matter interaction Hamiltonian in the RWA as 


H int (t) = £(t)V'+ £'(t)V 


(33) 


where V = Y2i= a /, /- l «.:|'0( c l i s now the de-excitation operator . 

The contribution to the signal comes from various Liouville space pathways starting 
with populations and coherences. Corresponding expressions can be obtained from the 
diagrams shown in the Fig. ((H) following standard diagrammatic rules jdf)| . The total signal 
is S[^l (lu; cl> 0 , 0) = (SpopO^; w 0 , C) + 5 r ( !; ( 1 ) j l 1 (co>; cj 0 , 0). It depends parametrically on the optical 
frequency cj 0 and the Rabi frequency Q = /i£ 0 /2h of the driving field. For a system initially 
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prepared in a population state (diagrams (al)-(a8) in Fig. 4)) we obtain (see Appendix B) 



( 34 ) 


and for initial coherences (diagrams (a9)-(a!2) in Fig. 4) we get 



+ (£(u> - u Q )n 

acl^cbG ca;ca (cj ^o) 4” & (^) l^bc/^cbQcb \ca ))p\ 


'ba 


ss 


(35) 


where the propagator Gki;mn( w) = ((kl\(ul — iC ) 1 \mn }} and / is the identity matrix in 


Liouville space. C is defined in Appendix B. In the expression for the signal the ojq depen¬ 
dence enters through the field as well as the retarded propagator. The presence of uj 0 makes 
the signal phase dependent. When cn 0 = fl = 0 we recover the equilibrium result (diagrams 
(al),(a3), (a5) and (a7)) and the signal reduces to 



(36) 


which solely depends on the power spectrum of the field and is independent of its phase. 

B. Simulations of the linear response of a strongly driven three-level system 

In the following numerical calculations we first solve for the steady state p ss using the 
master equation given in Appendix B (see Eq. fITTIi ) and then calculate the linear signal 
using Eq. ffTfll and Eq. (l35lh We use the linearly chirped Gaussian electric field S(uj) given 
in Eq. (Ell) . Fig 151(a) shows the equilibrium case when there is no driving i.e., ujq = G = 0- 
The signal then shows two peaks for the dipole transitions between the states |a) |c) 

and | b) — > |c). The peak at cu — uj ac is much stronger compared to the one at 0 Jt >c since 
Paa ~ Pcc Pbb ~ Pcc at equilibrium. Fig. [51(b) depicts the static situation with aq = 0 
but finite f2. This static coupling between states |a) and \b) renormalizes the energy values 
for these levels. These new energies are given by cu' a = |(w a + (Vb) — |(4Q 2 + ojl a )^ and 
bo' b = \{uj a + u>b) + |(4 Vt 2 + col a )^ with the gap between the two energy states is (4Q 2 + col a )^ 
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FIG. 5. Linear transmission signal S^\uj) for a strongly driven three level model system as a 
function of detected frequency co for two different cases, (a) - no driving = 0, c^o = 0, (b) - static 
coupling: c^o = 0, = 0.05 eV. Other parameters are: uo a = 0 eV, = 0.01 eV, oo c = 1.0 eV, 

ksT — 0.0259 eV, 7 ^ = 0.004 eV, y ca = 0.0001 eV, 7 ^ = 0.0002 eV. Parameters for electric field 
are Tq = 0.14 fs, O c = 0.5 eV, <fi" — 0. 
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Mb T 2 A a 6 2^ 

M a — ^Aafc + 

M a ~ §A a6 - 


Bare system system in strong field 


FIG. 6. (Color online) Frequency spectrum for the three level system with lower two levels driven 
by a strong field with frequency uq. Here A a b = mq — Mb a and = (40 2 + A 2 6 ) 2 . For details see 
Ref. 


a sir 

f- 


(see Fig. Elat ujq = 0). Here u>b a = Mb — M a . The signal shows two peaks at m = u c — uj' h and 
m = m c - J a . 

We next turn to the case of periodic driving with finite ujq. Fig. EKa) depicts the signal 
SH(M) generated by the system prepared initially in the population state given by a field 
with moderate strength G < cuo- In the presence of a monochromatic field the atomic 
wavefunction oscillates with four different frequencies js| (see Fig. El) : M a — |(A a b ±C') and 
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FIG. 7. Left column: Spop(w) - (a), - (b), (w) - (c) for moderate driving Q = 0.005 

eY, ujq = 0.01 eY. Right column: the same as the left but for the strong driving Q = 0.05 eV, 
wo = 0.01 eV. Other parameters are the same as in Fig. [5] 


u>b + |(A a b ± fl') where A a t> = w 0 — u>b a and A ac = w 0 — w ca are the detuning frequencies 
and Q' — (4fi 2 + A 2 J t To that end we observe four peaks in the signal which corresponds 
to following resonance frequencies w ca — |(A a & ± O'), and w c (, + |( A a j } ± O'). Similarly we 
obtain the signal for the system initially prepared in coherence. Fig. 01)) contains 

the same four peaks as in Fig. 0(a) but with the different intensity profile compared to 
the population contribution, due to different quantum pathways contributing to the signal. 
The overall intensity is suppressed, as coherence contribution is always weaker than the 
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(a) (b) 



Q (eV) Mi (eV) 


FIG. 8. (Color online) (a) absolute value of linear transmission signal |S f ( 1 )(a;)| for three level model 
system vs cu and Q at fixed c^o = 0.02 eV. (b) - the same signal but vs co and c^o at fixed Q — 0.1 
eV. Other parameters are the same as in Fig. [5] 

population. The total signal is depicted in Fig. 0(c) and is dominated by the 

population contribution. 

So far we discussed moderate driving Q' < cu 0 . In this case the gap between peaks 1 
and 2 and between peaks 3 and 4 is fl' whereas the gap between peak 1 and 3 and between 
2 and 4 is ujq (see Fig. ED- For strong driving fl' > ujq due to the level crossing of the 
energy spectrum the gap between peak 1 and 2 and between 3 and 4 becomes cn 0 whereas 
the gap between 1 and 3 and between 2 and 4 is Q!. The corresponding signals generated 
from the system initially prepared in the population state is shown in Fig. 0(d). The 
coherence contribution is shown in Fig. 0(e) and the total signal is plotted in Fig. 0 (f). 
Strong driving results in a higher frequency resolution in the spectra where all four peaks 
are well separated in frequency. Furthermore, the strong driving enhances the population 
contribution and suppresses the coherence contribution in the signal. 

The effect of level crossing and dependence of the signal with respect to ujq and Q for 
a broad range of parameters is shown in Fig. (El). Fig. El(a) depicts the absolute value of 
the signal |<S'( 1 )(o;)| vs the detected frequency u> and the Rabi frequency Q for fixed value of 
ujq = 0.02 eV. The signal clearly shows four peaks with fixed splitting equal to cn 0 between 
the nearest peaks (1 and 2) and (3 and 4). In the same time the gap between peak 1 and 3 
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FIG. 9. (Color online) (a) 2D contour plot of absolute value linear transmission signal |S'^ 1 - ) (a;)| vs 
detected frequency uj and chirp rate <fi" starting with initial population, (b) starting with coherence 
and (c) the total signal. The values of the parameters are coi,,, = 0.05 eV, uj ca = 1.0 eV, ujq = 0.05 
eV, = 0.15 eV, Co c = 0.3 eV. Other parameters are the same as in Fig. El 

and between 2 and 4 is O' which increases linearly with 0. Fig|8] (b) represents the signal 
vs uj and field frequency ujq for fixed value of Q = 0.1 eV. For small cuo the separation 
between the the peaks are not well resolved. However for higher uj 0 the separation between 
the nearest peaks increases linearly with ujq. 

In Fig. m we display the chirp rate dependence of the linear signal. If the system 
initially prepared in the population, the signal Spop shows oscillatory behavior for all four 
peaks with respect to the chirp rate due to the error function in the field 8 . The period 
of oscillations for the (1 and 2) pair of states with high energy is slightly shorter than for 
the (3 and 4) pair of states with low energy. This is due to the since central frequency of 
the pulse (ui c = 0.3 eV) which results in the faster oscillation of uj ca compared to uj c },■ This 
asymmetry is enhanced for the case when the system is prepared initially in coherence. Fig. 
Mb) depicts the 5'^, where in addition to error function, the nonlinear phase dependence 
comes from the oscillating exponent. However since, the total signal is dominated by 
the population contribution, the spectra in Fig. |9](c) is nearly indistinguishable from the 
population contribution alone. Similar population dominating behavior of the total signal 
was observed in Fig. 1. However unlike Fig.l where we did not address the details of the 
preparation in a driven system the initial conditions can be controlled by the driving field. 
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FIG. 10. (Color online) First row: density matrix elements vs c^o for ft = 0.01 eV: population 
Paa - (a), Pbb - (b), coherence Ke[p a b] - (c), and Im [p a b\- (d). Second row: 2D plots for the 
linear transmission signal |£ , ( 1 )(a;)| vs detected frequency uj and chirp rate (/)" starting with initial 
population - (e), starting with coherence - (f) and the total signal - (g) for ft — 0.01 eV and 
(jJq = 0.035 eV. Third row - the same as second row but for ft = 0.01 eV and c^o = 0.05 eV. Here 
ujba — 0.05 eV, uo ca = 1.0 eV. Other parameters are the same as in Fig. [5l 


In fact the steady state populations and coherences depend differently on the parameters of 
the driving field [8] such as ft and which makes it possible, in principle, to control relative 
contributions of populations and coherence and extract a pure coherence contribution. 

To demonstrate the effect of the control of the population/coherence contributions to 
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the signal we examine first the steady state density matrix elements as a function of D 
and uq. Fig. HHa) and (b) show the steady-state populations p aa and respectively as 
a function of Uq for the fixed value of the Rabi frequency fl. For a given parameters p aa 
( Pbb ) are symmetric functions of its argument which reaches its minimum (maximum) at the 
resonance frequency u 0 = 0.05 eV. Both populations are concentrated around the value of 
0.5 as the total populations of all states is always unity and p cc = 0 in the steady state. The 
same applies to the imaginary part of the coherence Im[p a &] as shown in Fig. ITUlfd) with 
the value changing from —0.01 to —0.07. On the other hand the real part of the coherence 
Re[p a b] is an asymmetric function of its argument ranging from 0.2 to —0.2 and reaches zero 
at uq = 0.05 with the maximum value around u 0 = 0.035 eV and u 0 = 0.055 eV. Therefore 
we have identified two parameter regimes when coherence contribution is substantial with 
the steady state value of 0.2 (Q = 0.01 eV, uq = 0.035 eV - parameter regime I) and regime 
when it can be neglected (fi = 0.01 eV, Uq = 0.05 eV - parameter regime II). Fig. [TTilf e) 
shows the oscillations of the population contribution with respect to nonlinear phase 0" using 
for parameter regime I. Due to weak driving Q = 0.01 not all four peaks depicted previously 
in Fig. [9] are well resolved. In fact there is one strong peak at u — 0.98 eV and two weak 
peaks at u = 0.94 eV and u = 1.0 eV. Coherence contribution to the signal is shown in Fig. 
[Turn . In this case all three peaks are manifested much stronger. Finally for the total signal 
shown in Fig. fTHlf g) we see that out of three peaks the strong peak at u = 0.98 eV and 
weak peak at u = 0.96 are suppressed whereas the weak peak at u; = 1.0 eV is enhanced 
compared to the pure population contribution shown in Fig. flOlf e). This is a manifestation 
of the strong coherence. In the parameter regime II the field frequency is u = 0.05 eV and 
the corresponding value of the steady state coherence drops by roughly a factor of 3. In 
this case the population contribution shown in Fig. fl(¥ hl dominates over the coherence 
contribution shown Fig. HDKi) and the total signal in Fig. fTTili ) is completely dominated by 
the population contribution. Therefore, a driven preparation allows to selectively manipulate 
the initial conditions of the system and separate the quantum pathways contributions from 
initial populations and coherences to the linear absorption signal. Note that similar analysis 
can be performed if we instead fix the value of u 0 and change the intensity of the driving 
field via fl. 
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V. CONCLUSIONS 


We have calculated the linear and non-linear frequency domain optical signals for systems 
prepared in a nonequilibrium state. The generalized n-tli order susceptibility x (n> then 
depends on n+1 independent frequency variables rather than n, which is the case for systems 
initially in equilibrium. This nonequilibrium state results in a nontrivial phase dependence 
of the electric field, already in the linear signal, even if the system is initially in a population 
state. This phase dependence is strongly enhanced by initial coherences. Furthermore 
we predict the new resonances in nonlinear signals that depend on the frequency difference 
between the initially prepared coherent superposition of molecular states. We then addressed 
a particular case of the initial preparation of the system. We investigated a three level 
system with two lower energy states strongly driven by a monochromatic field. In this 
case the driving field frequency cn 0 generates a phase dependent linear signal. Performing 
numerical simulations we show that depending on the strength of the driving field via its Rabi 
frequency, the relative contributions to the signals due to initial populations and coherences 
can be controlled and even a pure coherence contribution can be extracted. 


The new resonances and phase dependence demonstrated in the present work are closely 
related to quantum coherences. For instance, the atomic gas strongly driven with microwave 
radiation allows to alter detailed balance conditions and change the transmission properties 
of the optical pulses under certain conditions for electronic transitions and parameters of 


the driving and probe fields. These effects have 


:>een studied in cold alkali gases where ef¬ 


fects such as lasing without popu 


parency 


3l|, and slow light 32|, 


in semiconductor quantum dots 


ation inversion 


13-15], electromagnetically induced trans- 


33] were observed. These effects 


rave also been observed 


H 


35], and heterostructures 36|, 


37]. These fundamental 


effects are typically examined through simple linear transmission/absorption experiments 
well described by our formalism. Furthermore, the parameter regime for observation of co¬ 
herence effects is very similar to the one used in the simulations in Figs. 5-10. The power 
of the approach presented in this paper is not limited to linear experiments. The results of 
Eqs. (26)-(31) are described in terms of Liouville space electronic Green’s functions and thus 
allow to extend quantum coherence effects to higher order nonlinear optical measurements 
exploring interplay between coherence and many-body effects in atomic, molecular and solid 
state systems. These include electron and nuclear dynamics, nonadiabatic dynamics, chemi- 
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cal reactions and other many-body effects. We examined a series of measurements and their 
understanding in the context of spectroscopy of e.g. strongly driven systems. We demon¬ 
strated how electronic and vibrational structure and dynamics of molecules influenced by 
strong driving fields. The practical applications of our formalism in the context of quantum 
coherence effects discussed above has great potential in solid state technology. In particular, 
in the recent years quantum coherence in nitrogen vacancy (NV) defect centers in diamond 
show great promise for quantum information processing [2a], magnetometry 39-41] and 


electrometry [42 
ting constant |43|, 


. The recently measured temperature dependence of the zero-field split- 


44j indicates that it may also be used as an atomic temperature sensors. 


I . 


The broadband excitation of the NV centers in diamond with chirped pulses |45| allow the 
quantum control of the phases of single electron spins. The phase dependence of the optical 


signals shown by authors in Fig. 5 of Ref. 


45] yields similar features to our Fig. 7 of the 


present paper. The predicted resonances and phase control along with nonlinear optical 
measurements can be tested in these solid state systems. 
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APPENDIX A: FOUR WAVE MIXING SIGNALS 


In Fig. firm we draw all possible Liouville space pathways (for both initial populations 
and initial coherences) that contribute to FWM signal for the three level model system 


in Fig. (Oja. For diagrammatic rules we refer to Ref. 30]. The signal is calculated using 


Eq. m considering three incoming fields £{, i = 1,2,3 as monochromatic continuous waves 
and the probe field £4 as a spectrally broad pulses. The expression for the signal for these 
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FIG. 11. Ladder diagrams for three level system shown in Fig. (pjt). Here the states \i),\j),\k) 
representes two lower energy states | a) and \b) of three level system. 


pathways is given as 


S < £ > (u-,cj 1 ,uj2,uj 3 ) = -Re ^2 Pij 


i : j,k£a : b 


S* 4 (cu)£\Q£ x 


x 


l^jc^kc^kc^ic 


1 


X- 


CJ + CJ 2 —^3 — bJ c j-\~iT) 


{bJ ~ bJ\ +6^2 — ^3 — ^ij (cj — UJ c j -^rlTj) (bJ — CJ 3 — LU^j -\~iTj) 

+ (bJi bJ 3 ) (37) 


Sja (.d)-,u>i,u>2,u) 3 ) — -Re ^2 Pij 


h 


X- 


i,j,kEa,b 


£^(uj)£ 3 £*£ 1 


x 


l^kcl^cjl^ kd^ic 


1 


CJ+CJ 2 —^3 — bJ c j-\-if] 


(o;-cji+cj 2 -^3-^ + *??) (w-o; ci: +i?7)(w-a;3+zr]) 

+ (wi < > W 3 ) (38) 


S^(uj]UJ 1 ,UJ2,Ld 3 ) = -Re ^2 PH 

i,j,k£a,b 


X- 


£* 4 (cu)£ 3 £%£ 1 


x 


Pkcl^cjji'kcPic 


co—oji —Lu 3 —u)i c +irj 


(cu-u}i+u> 2 -co 3 -cUij + ir]) (cu-u> ck +iri)(cu-u>i-cuik+irj) 

+ (^i ^ ^ 3 ) (39) 


/o\ ^^ 

SaA <^2,^)3) =-Re 2 _^ Pij 


h 


z,j,fcGa,6 

1 


£l(u)£ i £*£ l 


X 


l^kcl^cj l^kc^ic 


X- 


CJ —CJi —bJ^ — bJic + ir] 


(cu-Ui+^-uJs-uJij + irj) (u!-u> cj +iri)(LO-u) 3 -\-iri) 

+ (t^i -f->- W 3 ) (30) 


The total signal is S®,, = S<? + Sf + Sf + Sg’. 
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APPENDIX B : MASTER EQUATION FOR A DRIVEN THREE LEVEL SYS¬ 
TEM AND THE LINEAR RESPONSE SIGNAL 


In this section we give the details about the master equation for the driven three level 
system described by the Hamiltonian given in Eq. (1321) . We follow the standard set of 
approximations such as weak system-bath coupling (second order), large reservoirs, wide- 


30, 


46] to derive the equations of motions for populations and 


band of the leads 8, l], 
coherences. 

We Define new set of variables (denoted by the symbol tilde) in the rotating frame 
Pab(t ) = e~ luot p ab (t), p ac (t ) = e~ luJot p ac (t ), p bc = p bc and p u = p ih i = a,b,c , to obtain the 
Bloch equation 


Paa i-^(Pba Pab) ('~i'oh '~]iac)Paa 4” ‘'tbaPbb 4” ^fcaPcc 

Pbb 11 ( f) ba Pab) (Pfba 4 “ 7 fec) Pbb 4 “ / ~fabPaa 4 “ / ~fcbPcc 

Pab ‘i^abPab 4 “ ^^-ipbb Paa) ^ (haft 4 “ Ota 4 “ 7 ac 4 “ O'tc )Pah 

Pac i^cicPeic ~\~ Tl/fo-. ^(0ac 4” 7ca 4” 0 ab 4“ '~fcb)Pac 

Pbc i(pJb ^c)Pbc 4“ i^Pac ^(bl>c 4“ 7c 6 4“ Ota 4“ 0V.:a) Pbc (41) 

where 7 ij,i,j = a,b,c are the decay rates from the state i to state j. Note that, the driving 
field couples the population of the lower two states with its coherence. The crucial advantage 
of working with the rotating frame is that the Liouville operator in this frame C becomes 
time-independent and the steady-state solution can be obtained. We write the above equa¬ 
tion in matrix form as p(t) = Cp{t) where p T = (p aai Pbb, p cc , Pab, Pba, Pac, Pea, Pbc, Pcb)- The 
transformation matrix between lab and rotating frame is diagonal 


U(t) = diag(l, 1,1, e iuot , e~ iuJo \ e iuJot , e~ iuJot , 1,1) 


(42) 


where diag represents the diagonal elements of the matrix. Since p{t) = U(t)p(t ) it implies 
that the propagators in the two frames are related via 

g{t,t') = u(t)g(t-t')u- l (t') ( 43 ) 

where Q(t,t') = —1 9(t — t , )e^ ^_^, ) and 1 9(t — t r ) TeJ)' drC< - T '> for t > t'. We first 

solve for the steady state in the rotating frame by demanding p ss = 0. Then in the lab frame 
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the solution for non-stationary state p nss at a particular time To (taken as the preparation 
time) is given as p nss (T, o) = U(ro)p ss which also mean that the coherence elements oscillate 
with the driving field frequency ujq even in the long-time limit. 

We probe the system, prepared in the state p nS s( T o ), with a weak probe field £{t) which 
allows transitions between the states |a) —» |c) and |6) —>■ | c). We impose the RWA and 
write the light-matter interaction Hamiltonian as 

H int (t) = £(t)V^+£*(t)V (44) 


where V = Yhi= ab Pi C \i)(c\. Given the initial state p nss (ro) we write the expression for the 
linear signal following Eq. flSD as 



POO 

rt 

ih£*(u) 

/ dte iut / 


J —oo 

J —oo 


coo, O) = — Im 

n _ 

(45) 

Using the relation between the propagators Q(t , t') and Q(t — t ') given in Eq. (T431) we obtain 


— — Im 

h 


£*(co) 


dte iujt f dT l £(T l ){V L U(t)g(t-T l )U-\T l )Vl_U(T l )) Ps 
J — oo 


(46) 


Note that in the last line the average is expressed with respect to p ss which is obtained by 
using the relation Q{j\ — To)p ss = —\pss- Performing Fourier transformation for the field 
and the propagator and writing down the matter correlation function explicitly by reading 
the diagrams for population and coherence in Fig. 0] we obtain Eq. (1341) and Eq. (J35J). 
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